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Abstract. Elliptic flow measurements as a function of pt of charged {tt^ and low-p^- 
protons) and strange (A and K^) particles from Pb+Au collisions at 158 AGeV/c 
are presented, together with measurements of and meson production. A 
mass ordering effect was observed. Scaling to the number of constituent quarks 
and transverse rapidity y^p^ scaling are presented. The results are compared with 
results from the NA49 and STAR experiments and with hydrodynamical calculations. 
For the first time in heavy-ion collisions, cj) mesons were reconstructed in the same 
experiment both in the K^K~ and in the e^e~ decay channels. The obtained 
transverse mass distributions of 4> mesons are compared with results from the NA49 
and NA50 experiments. The yield and the inverse slope parameter of the mesons 
were reconstructed from two independent analyses. Our results are compared with 
those from the NA49 and NA57 experiments. 
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1. Introduction 

Elliptic flow is described by the differential second Fourier coefficient of the azimuthal 
momentum distribution V2iT>) = (cos(20))x) [11121 El- The brackets denote averaging 
over many particles and events, and V represents a phase-space window in the {pT,y) 
plane in which f 2 is calculated. The azimuthal angle is measured with respect to the 
reaction plane defined by the impact parameter vector b and the beam direction. For 
non-central collisions (6 7^ 0), V2 is an important observable due to its sensitivity to the 
EoS, and through it to a possible phase transition to the QGP. Since we could identify 
protons via dE/dx only at low px-, the V2 of the A is important because this is a baryon as 
well. In comparison to the elliptic flow of pions and mesons the A flow can be used 
to check the mass ordering effect and for comparison to hydrodynamical predictions. 
Testing the differential flow measurements of different particle species against different 
scaling scenarios may yield additional information about the origin of flow. 

As strangeness enhancement has been suggested as a signature of the deconfined 
stage |3], understanding of the and meson production is important as here hidden 
and open strangeness are involved. The study of (j) yields in different decay channels is 
important in light of a possible modification of the mass, width and the branching 
ratios near the phase boundary. 

2. Experiment 

The CERES experiment consists of two radial Silicon Drift Detectors (SDD), two Ring 
Imaging CHerenkov (RICH) detectors and a radial drift Time Projection Chamber 
(TPC). The CERES spectrometer covers rj = 2.05 — 2.70 with full azimuthal acceptance. 
The two SDDs are located at 10 and 13 cm downstream of a segmented Au target. 
They were used for the tracking and vertex reconstruction. The purpose of the RICH 
detectors is electron identification. The new radial-drift TPC operated inside a magnetic 
field with a maximal radial component of 0.5 T providing a precise determination of the 
momentum. Charged particles emitted from the target are reconstructed by matching 
track segments in the SDD and in the TPC using a momentum-dependent matching 
window. A more detailed description of the CERES experiment can be found in 0. For 
the flow analysis, we used 30-10^ Pb+Au events at 158 AGeV/c collected in the year 
of 2000 data taking period. Of these, 91.2% were triggered on cr/ageo < 7%, and 8.3% 
events with a/ageo < 20%. The meson analysis in the kaon (dilepton) channel used 
24-10^ (18-10^) events taken with the most central trigger. 

3. Methods of strange particle reconstruction 

The A particles were reconstructed via the decay channel A — > p+7c~ with a BR = 63.9% 
and CT = 7.89 cm |5'^. Due to the late decay of the A particle, as candidates for A 
daughters, only those TPC tracks which have no match to a SDD track were chosen. 
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Partial particle identification (PID) was performed using dE/dx information from the 
TPC by applying a ±1.5o" (+lcr) window around the momentum dependent Bethe- 
Bloch value for pions (protons). On the pair level, a px dependent opening angle cut is 
applied, in addition to a cut in the Armenteros-Podalanski variables {qt < 0.125 GeV/c 
and < a < 0.65) to suppress Kg. With these cuts values for S/B ^ 0.04 and 
S/y/B ^ 500 were obtained 0. 

The Kg particles were reconstructed via the decay channel Kg ^ 7r+ + vr^ with 
a BR = 68.95% and cr = 2.68 cm [6 . Partial PID for 7r+ and tt~ was performed by 
applying a ±1.5cr window around the momentum dependent Bethe-Bloch energy loss 
value for pions. As the Kg particle comes from a primary vertex, a possibility to suppress 
fake track combinations is given by a cut (0.02 cm) on the radial distance between the 
point where the back extrapolated momentum vector of the Kg candidate intersects 
the X — y plane and the primary vertex. In addition, a cut of 1 cm on the z-position 
of the secondary vertex was applied. In this approach, the values of S/B ^ 0.92 and 
S/y/B ^ 500 were obtained |H1II|. 

In order to remove the effect of autocorrelations, tracks which were chosen as 
candidates for daughter particles were not used for the determination of the reaction 
plane orientation. In the case of A particle reconstruction, the combinatorial background 
was determined by ten random rotations of positive daughter tracks around the beam 
axis and constructing the invariant mass distribution, while in the case of Kg particle 
reconstruction, the mixed event technique was used. 

A (Kg) particles were reconstructed in y-pT-(p bins. We used the area under the 
peak, obtained by fitting the invariant mass distribution with a Gaussian, to measure 
the yield of A (Kg) in a given bin. Plotting the yield versus for different pt and 
y values one can construct the dN/^(^j^o-j/d(l) distribution. Fitting these distributions 
with a function c[l + 2^2 cos(2(/))], it is possible to extract the observed differential v'2 
values. The obtained f 2 coefficients were corrected for the reaction plane resolution via 
V2 = V2/ -y/2(cos[2($Q — $b)]) [3]. Here, $a and ${, denote the azimuthal orientations 
of reaction planes reconstructed from two random subevents. In the case of the vr^ 
elliptic flow analysis, subevents are formed from positive and negative pions separately. 
Using the method of subevents, correction factors were calculated for different centrality 
bins. In all 3 analyses (A, Kg and vr^) similar values were obtained. The corresponding 
resolution ranges from about 0.16 to 0.31, depending on the centrality. 

Due to the small statistics of strange particles, the differential elliptic flow analysis 
was performed for only two centrality classes. The huge statistics of tt^ allowed to 
perform the differential elliptic flow analysis in six centrality bins. As we used the 
combination of data taken with different triggers, the centrality is characterized by a 
weighted mean centrality (^^) calculated using the numbers of TPC tracks as statistical 
weights 0. 
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4. Results 



In Fig. n are shown the resulting px dependences of f 2 for three particle species. An 
increase of the elliptic flow magnitude vs pt for all three particle species is visible. In the 
case of A elliptic flow, the absolute systematic error Av2, estimated from two different 
ways of A reconstruction, is ^qq^^ for px < 1.6 GeV/c and ^qq2 Pt > 1-6 GeV/c 
which is small compared to the statistical errors. Particles are accepted as vr^ if their 
TPC dE/dx is within a ±1.5o" window around the nominal Bethe-Bloch value for pions. 
The HBT contribution to the vr^ elliptic flow is subtracted using the procedure described 
in jO]. Separately calculated elliptic flow of 7r+ and tt^ shows that the averaged difference 
between them is ~ 0.003 in both 77 and y, which can be attributed to the contamination 
of protons in 7r+ sample. Comparing results obtained from two independent analysis 
methods we concluded that the overall absolute systematic error in tt^ elliptic flow 
measurements is not bigger than 0.0036. 
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Figure 1. The A (left), (middle) and tt^ (right) elliptic flow vs transverse 
momentum in semicentral events. Hydrodynamical predictions are presented for two 
freeze-out temperatures: T/ = 120 MeV (sohd) and T/ = 160 MeV (dotted). 



The elliptic flow results are compared with the hydrodynamical calculations done by 
P. Huovinen based on |10t lllj. The calculation was done in 2+1 dimensions with initial 
conditions flxed via a fit to the pt spectra of negatively charged particles and protons in 
Pb+Pb collisions at 158 A GeV/c The underlying EoS assumes a first order phase 
transition to a QGP at a critical temperature of Tc = 165 MeV. The hydrodynamical 
predictions were calculated with 2 freeze-out temperatures, Tj = 120 MeV and Tj = 
160 MeV. The model prediction with the lower freeze-out temperature of Tj = 120 MeV 
overpredicts the data, while rather good agreement can be achieved with a higher freeze- 
out temperature of Ty = 160 MeV (this is however not the preferred value considering 
the proton pt spectra). 

A comparison of the CERES data to results from NA49 JHl at the same energy 
(^5^ =17 GeV) and to STAR results dH at = 200 GeV is shown in Fig. H 

The NA49 and CERES data are in very good agreement. After rescaling the STAR 
results to the centrality used in the CERES experiment, the V2 values measured at 
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Figure 2. Comparison of A (left) and Kg (middle) elliptic flow measured by CERES, 
STAR and NA49. Comparison between the elliptic flow magnitude of the tt^, low-pT 
protons, A, and Kg in semicentral events (right). 



RHIC are 15 — 20% higher due to the higher beam energy. In Fig. |2l (right), the elhptic 
flow magnitude of the vr^, Kg, low momentum protons, and A measured by CERES 
are compared. A mass ordering effect is observed. At small pt, up to ~ 1.5 GeV/c, 
f2(A) < V2{Kg) < V2{7!'^). In the region of high pt, above ~ 2 GeV/c, the tendency is 
the opposite. As proton and A hyperon have similar masses and 3 valence quarks each, 
the f 2 of low momentum identified protons is considered as a natural continuation of A 
V2{pt) dependence in the region of small pt- The indication of a possible undershoot to 
negative values is tantalizing but not significant in view of the statistical errors. 
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Figure 3. Comparison between the elliptic flow magnitude of tt*, low-pT protons. A, 
and Kg scaled to the number of the constituent quarks (left) and to the y^^ variable 
(right). 



Fig. El (left) shows the scaled elliptic flow magnitude V2/nq for 7?=*=, K^, low-p^ 
protons and A plotted against pr/'nq in semicentral events. Here, Uq denotes the number 
of the constituent quarks. There is an indication that high p^ particles {pt > 1-5 GeV/c) 
show scaling behavior. A similar behavior is observed by the STAR experiment at RHIC 
|14j . This is consistent with the coalescence mechanism where co-moving quarks with 
high Pt form hadrons. In this case scaling to the number of the constituent quarks 
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shows the original momentum space azimuthal anisotropy formed at the early stage of 
the collision. 

Within the Buda-Lund model of hydrodynamics ^H], a scaling of elliptic flow of 
different particle species has been suggested jXTHinj when instead transverse momentum 
the transverse rapidity is used. We use their scaling variable yjf jlH] and show, in Fig. El 
(right), the results for tt^. Kg, low-pT protons and A in semicentral events. Within 
statistical errors a reasonable scaling is observed for all particles. This may indicate a 
hydrodynamic behavior of matter created in central heavy-ion collisions at the highest 
SPS energy. 
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Figure 4. Transverse momentum and rapidity Kg spectra from the Kg analysis 
performed without PID and without secondary vertex reconstruction |19| . 

Two independent analyses of the Kg spectra were done using the CERES 
data IS]. The first one, performed without PID and without secondary vertex 
reconstruction, is based on TPC information only [19 . A cut in the Armenteros- 
Podalanski plane was used in order to suppress A contamination. The pt and y spectra 
are shown in Fig. |3] An alternative approach of the Kg reconstruction was performed 
without PID but with secondary vertex reconstruction jH] which is already described in 
Section El In both analyses, Pb+Au events taken with the most central trigger were 
used. The Kg transverse momenta spectrum obtained with this analysis [2, is shown 
in Fig. El (left). The invariant multiplicity was fitted with an exponential fall-off with 
transverse mass mj. The yields and the inverse slope parameter T of the pt spectra 
from the two analyses are in good agreement. A comparison with results from other 
experiments is shown in Fig. (right). In order to match the centrality of the CERES 
experiment, results from the NA49 fK)\ E] and NA57 experiments are slightly 
rescaled. A rather good agreement between the NA49 analysis of charged kaons and 
the CERES Kg results in shape and yield was found. The difference in the yield is only 
5%. The rapidity distribution of Kg observed by NA49 shows a similar shape as the one 
from CERES (represented with the blue dotted line in Fig. El (right)) and a relatively 
good agreement in the yield. Within the CERES acceptance the results agree with the 
NA57 data, although the NA57 fit does not. 
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Figure 5. Left: Transverse momentum Kg spectrum from the Kg analysis performed 
with secondary vertex reconstruction j8j. Right: A comparison between CERES results 
(red circles |S] and blue squares published (open triangles) j^Ol and preliminary 

(open crosses ) 121] NA49 resuhs and NA57 data (green diamonds) g^. The black 
dotted line represents a fit to the charged kaon yield measured by the NA49, while 
blue (green) dotted line corresponds to a fit to the Kg yield measured by the CERES 
(NA57). 

The CERES experiment enabled for the first time at SPS to study simultaneously 
the leptonic and charged kaon decay modes of the (p meson, which may shed light 
onto the puzzle [23] • In order to obtain the Pt spectrum of mesons, the invariant 
mass distributions of K^K~ pairs were constructed. The corresponding distributions 




Figure 6. Left: K~^K invariant mass spectrum after background subtraction in 
1.5 GeV/c < < 1.75 GeV/c and 2.2 < < 2.4. Right: e+e" invariant mass 
spectrum compared to the hadron decay cocktail (solid line) and to a model calculation 
assuming the dilepton yield form the QGP phase and an in medium spread p (dashed 
line). 

of the combinatorial background were calculated using the mixed-event technique. 
An example is shown in Fig. IHl (left). To study mesons in the dilepton {e~^e~) 
decay mode, electrons are identified using the RICH detectors and the TPC dE/dx. 
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The main difficulties of reconstructing the meson in the dilepton channel are the 
low branching ratio and huge combinatorial background. Details of how to reduce 
the combinatorial background are explained in El ES]- The e~^e~ invariant-mass 
spectrum, corrected for the efficiency and normalized to the number of charged particles 
in the acceptance is shown in Fig. IHl (right). In the same figure are shown the 
expectations from the hadron decay cocktail [26j, as well as a model calculation where 
the cocktail p contribution is replaced by an explicit in-medium modification combined 
with continuous vrvr annihilation |2ZI. The later accounts very well for the data. The 
inverse slope parameter of T =273±9(stat)±10(syst) MeV and a rapidity density dN/dy 
of 2.05±0.14(stat)±0.25(syst) in the K+K~ mode and T =306±82(stat)±40(syst) MeV 
and dN/dy =2.04±0.49(stat)±0.32(syst) in the dilepton mode are in good agreement 
within errors. The data do not support a possible enhancement of the yield in the 
dilepton over the hadronic channel by a factor larger than 1.6 at the 95% CL. 




p, (GeV/c) m,-ms (GeWlc) 

Figure 7. Left: Acceptance and efficiency corrected px spectrum of (p measured in 
the K^K^ (open cycles) and e+e~ (closed circles) decay modes. Right: Scaled ttlt 
distribution of (p mesons reconstructed in the K'^K~ (triangles) and e+e" (circles) 
decay channels compared to the results from NA49 (squares) and NA50 (diamonds). 

The pt dependence of the meson yield measured in the K^K~ and e"''e~ decay 
channels, corrected for the acceptance and efficiency, is shown in Fig. [3 (left). The results 
are in very good agreement. After accounting for the slightly different measurement 
conditions, a comparison between CERES results and the existing Pb+Pb systematics 
j2H] is shown in Fig. [7| (right). The CERES results are in good agreement with the 
results from NA49 measured in the kaon channel. On the other hand, CERES data in 
the K^K~ channel do not agree with NA50 results. 
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